Superdiffusive nonequilibrium motion of an impurity in a Fermi sea 
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We treat the nonequilibrium motion of a single impurity atom in a low-temperature single-species 
Fermi sea, interacting via a contact interaction. In the nonequilibrium regime, the impurity does 
a superdiffusive geometric random walk where the typical distance traveled grows with time as 
~ f or ^-dimensional system with d > 2. For nonzero temperature T, this crosses over to 

diffusive motion at long times with diffusivity D ~ j 1 -^ -1 )/ 2 _ These results apply also to a nonzero 
concentration of impurity atoms as long as they remain dilute and nondegenerate. 
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In condensed matter physics, the behavior of a single 
impurity immersed in a sea of majority particles has been 
one of the simplest many-body problems and attracted 
much attention (see, e.g. [1, 2]). Especially, transport of 
an impurity in Bosonic [3] and Fcrmionic [4, 5] quantum 
liquids has been studied since a few decades ago. Re- 
cent realizations of ultra-cold polarized two-component 
atomic Fermi gases [6-12], with their remarkable con- 
trollability over parameters such as interaction strength 
and individual populations, have made it possible to 
directly access this type of quantum many-body sys- 
tem. Meanwhile, there have been a number of theoreti- 
cal works on the single impurity problem in an ultracold 
Fermi gas [13-17], investigating the equilibrium and near- 
equilibrium properties. Quite naturally, transport phe- 
nomena have also been studied [18] . More recently, trans- 
port experiments of an ultracold mass-balanced polarized 
Fermi gas [19] and an ultracold mass-imbalanced mixture 
with unequal populations [20] opened new opportunities 
to directly investigate nonequilibrium and dynamic prop- 
erties of population-imbalanced Fermi systems. 

In this paper, we discuss nonequilibrium transport in 
the high polarization and low temperature regime, with- 
out restricting ourselves within small deviations from 
equilibrium. As a limiting case of high polarization 
and low temperature, we first consider a single minority 
atom moving in a zero-temperature Fermi sea of majority 
atoms. We find that the impurity does a superdiffusive 
geometric random walk, where the time between colli- 
sions grows in proportion to time, and the impurity loses 
a fraction of order one of its energy in each collision. In 
d dimensions (d > 2) the typical distance traveled grows 
as ^/(rf+i). 

As is conventional, we call the majority species "up" 
f and the minority (impurity) species "down" J,- Note 
that in the regimes we study in this paper, the statistics 
of the minority atoms do not enter, so they may equally 
well be bosons or fermions. Assume the minority atom is 
initially near the origin in real space, with some proba- 
bility distribution of its momentum with <C /c^, 
where kp^ is the Fermi momentum of the majority Fermi 
sea. The minority atom is "dressed" either as a polaron 

or as a molecule, with effective mass m* and thus energy 

h 2 Q 2 

E^ 2 ; we choose the rest energy of the dressed mi- 



nority atom as its zero of energy. We assume that E± is 
low enough so that no internal excitations of the polaron 
or molecule are possible. Note that Fermi polarons and 
bosonic molecules in this high-polarization limit were ex- 
perimentally studied in Ref. [21]. 

We will treat the majority atoms as noninteracting, 
although the results we obtain appear to remain qual- 
itatively correct even if the majority atoms do weakly 
interact with one another, as long as the majority atoms 
remain a Fermi liquid. In the latter case, the minor- 
ity atom produces quasiparticles and quasiholes when it 
scatters from the majority Fermi liquid. We will first 
treat the case of a three-dimensional Fermi gas, return- 
ing to general d later. 

This dressed minority atom can scatter from a ma- 
jority atom with momentum slightly below kp-f and 
thus emit a majority particle-hole pair, which lowers the 
energy of the minority atom. Here we are considering 
such scattering events that occur "on-shell" ; the virtual 
particle-hole pairs that dress the polaron or molecule are 
assumed to have already been included and have renor- 
malized the rest energy and the effective mass of the 
dressed impurity atom. After the scattering event, the 
emitted majority particle and hole move away ballisti- 
cally at the majority Fermi velocity. In such a scattering 
event the dressed minority atom typically loses a frac- 
tion of order one of its energy. The rate of scattering is 
proportional to E? for d = 3, as we derive below. As a 
result the typical energy evolves with time as 



which results in E±(t) ~ t x / 2 and typical speed v ~ 
Ql ~ < -1 / 4 . Thus the minority atom does a superdif- 
fusive geometric random walk, where the typical time r 
between scattering events grows as r ~ t and the mean 
free path I grows as t ~ vt ~ t 3 / 4 . The last step of this 
walk typically gives the dominant contribution to the to- 
tal distance traveled, so the latter also grows as ~ i 3 / 4 . 
The number of majority particles and holes produced and 
the number of steps in this unusual random walk grow 
only as ~ logt. 

In order to understand this behavior more quantita- 
tively, we consider the time-dependence of the momen- 
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turn distribution /(Qj,,t) of the dressed minority atom. 
f(Ql,t) has a dimension of (length)^ and is normalized 
to unity when integrated over wavevectors Q^. We es- 
sentially treat the minority atom classically. Although its 
initial state may be a quantum-coherent wavepacket and 
thus of low entropy, as it produces particle-hole pairs it 
becomes more and more highly entangled with these ex- 
citations that it has produced in the majority Fermi sea. 
As a result, at long times its reduced density matrix is 
a mixed state, any initial coherence is transferred to the 
Fermi sea, and the volume occupied by the minority atom 
in its position-momentum configuration space grows as 
~ (iQl) 3 ~ t 3 / 2 , so its entropy grows as ~ logt. 

A single minority atom in an infinite system docs not 
disturb the distribution of the majority atoms, which are 
thus assumed to remain in a Fermi-Dirac distribution 
at all times. We now consider both the case of a zero- 
temperature majority Fermi well as low nonzero 
temperature T <C Tp-\ where the superdiffusivc behavior 
discussed above crosses over to diffusion at long time as 
the minority atom thermally equilibrates with the Fermi 
sea. The scattering process is a majority atom with mo- 
mentum q-f scattering from the dressed minority atom 
with momentum in to solid angle dCl in the rest frame 
of the Fermi sea, resulting in momenta and re- 
spectively, and the reverse of this process. The resulting 
time-dependence of the minority atom's momentum dis- 
tribution is 




dt 



d 3 



da 



v r — dn{n(e k ^f(Ql,t)(l-n(e qr )) 



FIG. 1: (Color online) Momentum space configuration of the 
four momenta described in the text. The lengths of vectors 
are not drawn to scale. 



the small energy transfer, \q\ ~ |A;| ~ ■ Using these 

axes to set up spherical polar angles, points along 
angles 9, cf>, where 9 is the angle from the z axis. With 
these coordinates, to the precision needed in this low en- 
ergy limit, the momentum transfer is (Q^ — QjJ = irjkp^. 
Fig.l schematically shows the configuration of the four 
momenta and the axes. We scale the (small) momentum 
transfer, defining 7 so 7Q = rjkp^. We also set the zero 
of energy for the majority atoms at their Fermi energy, 
and scale the (small) distances q and k from the Fermi 

surface, defining p, so that t q = h k ^ q = ^3rV- With 
all these definitions and scalings, we obtain 



(2tt) 3 ' dfl 

-«M/(Qi,<)(i-«K))] 



(2) 



d/(Qi,*) h Q* da 



d(f>d(cos9)dp-fd'y x 



dt (27r) 2 rn* dtt , 

[n(e kt )f(Q[, t)(l - n(e qf )) - n(e^)f(Q h t)(l - n( Cfct ))] 

(3) 



with e fct 



h 2 Q 2 



^(27% cos6> - 7 2 ). Now the inte- 



where v r is the speed of the relative motion and n(efe) is 

fc2 . 2 

the Fermi-Dirac distribution with e k = the energy 
of a majority atom. 

In the low energy <C kp^ and low tempera- 
ture ksT <C £f\ limits, the scattering is all at rela- 
tive momenta near kp^ with small momentum transfer 
|Qr — Q4.I "C k F ^. In these limits, v r ps and are 
essentially constant and may be taken outside of the in- 
tegral. The differential cross-section is dependent on the 
internal structure of the polaron or molecule, for which 
there are only approximate results, so we will just take 
it as an input parameter to our results. Near unitarity 
k F ^\a\ > 1 and |g 

k F ^\a\ < 1 and g§ 
scattering length. 

To scale the dynamics, we assume that distribution of 
Qj, has a characteristic momentum Q. Near equilibrium 

~ kgT, while for a "hot", nonequilibrium state 

Y% ^> ksT. In either case, the energy transfer in the 

scattering is of order \ <C e_Ff. Let's put the x axis 
along q^, so q-j- = ~x.(kp^ + q). The momentum transfer 
is small, so the angle rj between q^ and k^ is small 77 ~ 
1^ <SC 1; we define 77 to be nonnegative. Put the z axis 
so k^ = (x cos 77 + irj)(kp^ + k). Due to ksT <C ep\ and At zero temperature, the integral over p is readily done, 



, while away from unitarity 
2 , where a is the bare s-wave 



C 9T 1 2m* / Q 

gration measures are dimensionless, so the scaling of the 
scattering rate ~ Q 4 is explicit. Note that the result does 
not depend directly on the majority density, although 
near unitarity there is a dependence via the differential 
cross-section. 

If the distribution is isotropic in momentum space, 
f{Q\i t)> then it remains isotropic. If it starts anisotropic, 
the dynamics brings it asymptotically to isotropy at long 
times. Assuming isotropy, the integral over <f> can be 
done, and the integrals over 9 and 7 can be exchanged 
for an integral over Q'^ by multiplying J dQ^8(Q'^ — 

Ql + 7 2 <3 2 - 2Q;Q 7 cos 9), giving 
df(Qj,t) „ HQ 2 da f Q\ 

[n(e kt )f(Q[, t)(l - n(e qf )) - n(e qt )f(Q i ,t)(l - n(e kt ))}. 

(4) 
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and for the loss term also the integral over Q'^, giving 
df(Q^t) h da 



dt 



2irm* dQ 



( / dQ[Q[{Q^ - Ql)f(Ql t)) - -QiHQi, t) 

(5) 

Thus we see that at zero temperature, the total rate of 
scattering of a minority atom with momentum by 
producing a majority particle-hole pair and going to any 
lower energy is j^^Qf- 

This zero-temperature dynamics has a superdiffusive 
long-time scaling form in terms of a scaled momentum 



s = (At) 1 ^ 
and a scaling function g(s), with 

f(Q x ,t) = {Atf' A g{s) 



(6) 



(7) 



With A = „ . -g^ . the scaling function satisfies the 

27rm* d\l ' & 

integro-differential equation: 
3 , . s dq(s) 



-^s i g(s)+J^d S 's'(s' 2 -s 2 )g(s') 



(8) 



The resulting g(s) appears to be smooth and of order one 
at small s. Its large-s asymptotics are: 



s dg(s) 
4 ds 
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s 4 g(s) => g(s) 



exp 
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(9) 



A small "cloud" of minority atoms in a zero tempera- 
ture Fermi sea will also spread in this same superdiffusive 
fashion if it is effectively nondegenerate and noninteract- 
ing. In order for the minority atoms to be effectively 
nonintcracting, the majority particles and holes produced 
must leave the cloud without additional scattering, so the 
majority mean free path within the cloud must be larger 
than the cloud. If the minority cloud is initially too dense 
and/or too "hot" for this to be true, it will rapidly ex- 
pand and cool and cross over to this superdiffusive regime 
of behavior at later times. 

If the majority Fermi sea is not at zero tempera- 
ture, the superdiffusive behavior will continue until the 
minority atom cools down to near equilibrium, where 
^-S- ~ k B T. After that it will move diffusively, with 

typical speed v ~ y and scattering rate 1/rj, ~ 

AQ 4 ~ — ^f T ^ g^- The resulting equilibrium spin dif- 
fusivity is D s ~ v 2 t^ ~ fi 3 /(m* 2 k B T^). In previous 
work, the same temperature dependence in terms of the 
equilibrium spin-drag relaxation rate was obtained for a 
mass-imbalanced system [22]. Again, if there is instead 
a "cloud" of minority atoms, this will be their behavior 
when they are dilute enough to be nondegenerate. 



It should be noted that the "on-shell" condition dis- 
cussed in the beginning of the present paper is crucial in 
the analysis. As pointed out in Ref. [23] in the context of 
heavy ion transport in He 3 , if the definite effective mass 
assumption is not satisfied, the temperature dependence 
of transport coefficients can be significantly modified. We 
can understand this from the observation that the mass 
dependence of the scattering time is ~ m/m* . There- 
fore, when m/m* <C 1, the lifetime of the impurity can 
be small and thus its energy uncertainty becomes large 
and the assumption of clastic scattering breaks down. As 
a result, our analysis may not be applicable if the mass 
of the minority atom is much greater than that of the 
majority atom. 

Next we discuss a general o?-dimensional Fermi gas for 
d > 2. The argument is parallel to that given above 
for d = 3. The time evolution of the minority mo- 
mentum distribution /d(Q^,£) is the same as Eq.(2) ex- 
cept for the changes in dimension: d 3 q^- — > d d q^- and 
(27r) 3 — > (27r) rf . The solid angle d£l now possesses d- 
dimensional hyperspherical form. For <C fcpt the 
scattering cross section and relative speed remain essen- 
tially momentum independent and may be taken out- 
side of the integral. We use the same configuration 
of momenta (Fig.l) and dimensionless parameters p, m 
and 7 introduced above. Once we fix the direction of 
q^, we are free to perform a rigid body rotation to 
the system of four momenta in (d — 2) angular direc- 
tions. Thus, dQ = Sd-2(si~nri) d ~ 2 dri rj Sd-2V d ~ 2 dr] = 
Sd~2{Q/kFt) d ~ 1 J d ~ 2 d'~f (?? is small). 5^-2 is the surface 
area of a (d — 2)-sphere, given by Sd — 2n d / 2 /T(d/2), 
where T(x) is the gamma function. The integration mea- 
sure d d q^ is scaled as {m/m*)k <1 ^ 2 Q 2 dpdVt d , where dVl d 
is the solid angle element of on the hypersphere. We 
set the domain of both rj and 9 to be [0, 7r] for all d > 2. 
Although rj and 9 may be defined from to 2tt for d = 2, 
the reflection symmetry ensures us that we may integrate 
only from to ir and the factor of 2 can be absorbed in 
Sq = 2. With all these scalings, we obtain 



dfd 
dt 



do 

dn 



Sd- 

(27T) 



d+l 



j d - 2 d-fdn d d P [.. 



( 10 ) 

where we abbreviated the repeating phase fac- 
tor [n(e kf )f d (Q' v t)(l - n(e qf )) - n^)/^, t)(l - 

n ( e k t ))] — [■••]• {do/dQ)d is the generalized differential 
scattering "cross section" in d-dimensions. Now we can 
explicitly see the scaling of the scattering rate is ~ Q d+1 
in a nonequilibrium regime and thus - T^ d+1 ^ 2 near 
thermal equilibrium. The typical distance traveled / is 
~ £<V(d+i) an£ j -j- ne vo i ume m position-momentum config- 
uration space that the minority atom occupies increases 
as (lQi) d ~ . The growth rate of the number 

of scattering steps remains ~ logi. 

Assuming isotropy, we can replace dfld in Eq.(10) with 
Sd-2 (sin 9) d ~ 2 d9. As we did for d = 3, we multiply 

/ dQ[d(Q[ - Jof+^Q 2 - 2Q i Q 1 cos9). Then, for an 
odd d we do the 9 integral first, while for an even d we 
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do the 7 integral first. This gives 

df d (Qi,t) = h (da\ (Sd-2) 2 
dt m* \dfl J d (27r) d 

i3n' , 



dp[- 



Q-'Q'.dQ'.dj . n9 t , 

4- 



Q 2 Q[dQ[d9 
Qf-Q 2 sm 2 i 



7o d - 2 (sin0) 



d— 2 



(11) 

(12) 



where Eq.(ll) is for an odd d and Eq.(12) is for an even 
d. 70 is the value of ^(9) that satisfies the delta function 
(it has two roots and which root to use depends on the 
magnitudes of Qi, Q\, and 9). For an odd d, the integra- 
tion range of 7 is from \Q± — QjJ/Q to (Qi + Q'±)/Q and 
the integrand is a polynomial. Therefore, the 7 integral 
is elementary, although a formula for a general odd d is 
fairly lengthy so we do not present here explicitly. On the 
other hand, the 9 integral for an even d is not elementary 
so it may have to be evaluated numerically. 



We introduce A d 



h (do\ (Sd-2) 



which has di- 



mensions of (length( d+1 )/time). When d = 3, A d = A = 
27rm* as above. At zero temperature, we can do the 
integral / dp[...] to obtain 



dt 



(13) 



A,, 



dQ'MQ' i ,Q±)f d (Q , l ,t) - C d Q d + l f d {Q^t) 



The function hd(Q'p Q\) and the constant C d should be 
determined by Eq.(ll) or Eq.(12) and the p integral. For 
an odd d, C d can be found exactly, while for an even d, 
C d could be found numerically. As we saw earlier, for 
d = 3, C 3 = 2/15 and h 3 (Q[, Q x ) = Q[(Qf - Q 2 ). 

Before we scale the long time dynamics of Eq.(13), 
let us estimate the spin diffusivity in d dimensions. 
Near thermal equilibrium, v 2 ~ (hQ/m*) 2 ~ (fcsT)/m* 
and the scattering time t is ~ 1/ (A d Q d+1 ). Thus, 
the spin diffusivity in d dimensions is D s ~ v 2 t ~ 



m*( d + 1 ^ 2 (k B Ty d - 1 )/ 2 (d<T/dn) d 



when the minority gas is 
nondegenerate and the majority gas is an ideal Fermi 
sea or a Fermi liquid. 

It is straightforward to obtain a scaling form once we 
introduce the following scaling of the momentum and the 
distribution: 



Sd = (Adty/^Q, 

fd(Q h t) = (Adt) d l^g d { Sd ) 
The scaled integro-differential equation is, 



(14) 
(15) 



d , . s d dg d (s d ) 

9d{Sd) + 77377 j 

a + 1 dsd 



V 



ds' d g d (s' d )h(s' d ,s d ) 



C d s d d +1 g d (s d ) 



(16) 



We see that the asymptotic form of g d (s d ) at large s d is 
~ exp[-C d sf 1 ]. 

For d = 2, which is the other experimentally accessible 
case: A 2 = ^ (%) 2 , C 2 - 0.45 and h 2 (Q[,Qi) = 

— ?J — The spin diffusivity near thermal 



equilibrium in d = 2 is D s ~ TO , 3/2(fcBT)1/2(dg/rfo)2 ■ 

So far, our calculation relied heavily upon the phase 
space integral. In d = 1, however, the reduced phase 
space greatly restricts the scattering process. The scat- 
tering of an impurity from a T = Fermi sea for 
Qi *C via creating a single particle- hole pair is for- 
bidden in d = 1 by energy and momentum conservation 
unless m* <C m, which is not true in our case. Therefore, 
a low-momentum moving impurity in d = 1 is stable at 
zero temperature against scattering at this order. This is 
why the results of the present paper only apply for d > 2. 
The dynamics of a heavy impurity in a one-dimensional 
Luttinger liquid due to other higher-order scattering pro- 
cesses is discussed in Ref. [24]. 

Finally, we consider an experimental procedure to pos- 
sibly observe the expected superdiffusive behavior. Ini- 
tially, the majority and impurity atoms are trapped by 
species-selective potentials and cooled to very low T. The 
impurity atoms are tightly trapped at the center of a ma- 
jority cloud, and thus they are initially degenerate, but 
with a lower Fermi momentum than the majority atoms. 
(Here we discuss the case where the impurity atoms are 
fermions; if the impurity atoms are bosons, the specifics 
of how to set up an appropriate initial condition will be 
different). Then, release the impurity atoms but not the 
majority atoms and observe the subsequent expansion of 
the minority cloud using in situ imaging. At first, the 
impurity cloud should expand ballistically since scatter- 
ing is highly restricted by both majority and minority 
Pauli blocking. Quickly the impurity cloud expands and 
becomes nondegenerate and superheated relative to the 
majority gas. At this point the impurities start doing 
the superdiffusive motion discussed in this paper, which 
should be reflected in a superdiffusive expansion of the 
minority cloud. This superdiffusive motion persists until 
the impurity atoms cool to near the temperature of the 
majority Fermi gas. Then the impurity motion crosses 
over to standard diffusion. Recently, the expansion dy- 
namics of initially localized impurities in a 1-dimensional 
Bose gas has been realized and imaged [25]. 

In conclusion, we have studied the nonequilibrium 
and near-equilibrium motion of nondegenerate impurity 
atoms in a low-temperature Fermi sea. In the nonequi- 
librium regime we analyzed the unusual superdiffusive 
random geometric walk performed by the impurity atom 
for systems with dimensionality d > 2. At nonzero 
temperature this crosses over to standard diffusion at 
long times. We find that the equilibrium spin diffusiv- 
ity D s ~ 7 1 -( d ^ 1 )/ 2 in this regime where the majority 
atoms are degenerate, while the minority atoms are not. 
This temperature dependence of the spin diffusivity is 
something that may be measured soon in experiments 
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on 2- and 3-dimensional Fermi gases. The superdiffusivc 
nonequilibrium behavior will be more of a challenge to 
explore experimentally, since it requires a hierarchy of 

h 2 Q 2 

three energy scales: ksT <C -g— r <C ksTp^, and thus a 
very cold majority Fermi gas. 
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